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PREFACE 


The  potential  benefits  of  a  large  space  structure  (LSS)  for 
communications,  power  generation,  surveillance,  astronomy, 
and  space  defense  are  unlimited.  Consequently,  the  problem 
of  controlling  an  LSS  while  it  orbits  the  earth  is  receiving 
increased  attention  in  both  the  military  and  civilian 
communities.  The  problem  is  greatly  complicated  by  the 
extreme  size  and  flexibility  of  the  LSS  and  the  difficulties 
in  testing  an  LSS  control  system  in  the  1-g  environment  of 
the  earth's  surface.  As  a  result,  most  of  the  LSS  control 
work  to  date  is  theoretical  and  based  upon  computer 
simulation.  This  paper  is  no  exception.  It  proposes  an 
algorithm  which  aids  the  controls  engineer  in  specifying  a 
configuration  for  sensors  (motion  or  rate  detectors)  and 
actuators  (force  or  torque  generators)  to  control  an  LSS 
modeled  as  a  linear  stochastic  system  (i . e . ,  a  computer 
model)  . 


Subject  to  clearance,  this  manuscript  will  be  submitted  to 
the  American  Institute  for  Aeronautics  and  Astronautics' 
(AIAA)  Journal  of  Guidance.  Control,  and  Dynamics  for 


consideration.  As  a  result,  the  paper  is  written  at  a 
technical  level  assuming  the  reader  has  graduate  exposure  to 
structures,  dynamics,  and  optimal  stochastic  control. 
Additionally,  the  paper  is  tailored  to  the  AIAA’ s  formatting 
requirements  which  differ  from  thd  ACSC  requirements.  These 
differences  are  most  noticeable  in  line  spacing,  margins, 
table  and  figure  format,  heading  format,  reference  citation, 
page  number  location,  and  paragraph  indentation.  Also,  EDC 
has  authorized  substituting  the  AIAA  required  abstract  for 
the  Executive  Summary. 
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ABSTRACT 
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This  paper  presents  an  algorithm  which  aids  the  controls 
engineer  in  specifying  a  sensor  and  actuator  configuration 
for  regulation  of  large  scale,  linear,  stochastic  systems 
such  as  a  large  space  structure  (LSS)  model.  The  algorithm 
uses  a  linear  quadratic  gaussian  CLQ6)  controller,  an 
efficient  weight  selection  technique  based  upon  successive 
approximation,  and  a  measure  of  sensor  and  actuator 
effectiveness  to  specify  a  final  sensor  and  actuator 
configuration.  This  configuration  enables  the  closed-loop 
system  to  meet  output  specifications  with  minimal  input 
power.  The  algorithm  involves  no  complex  gradient 
calculations  and  is  numerically  tractable  for  large  linear 
models  as  demonstrated  by  the  solar  optical  telescope 
example  in  this  paper.  Additionally,  the  algorithm  provides 
the  controls  engineer  information  on  the  important  design 
issues  of  actuator  sizing,  reliability,  redundancy,  and 
optimal  number. 
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1.0  INTRODUCTION 


The  advent  of  the  Space  Shuttle  makes  the  large  space 
structure  (LSS)  an  imminent  reality.  These  future  space 
structures  will  be  measured  in  kilometers  and,  of  necessity, 
will  be  lightweight  and  highly  flexible  (light  damping) . 
Standard  LSS  missions  will  include  power  generation, 
surveillance,  astronomy,  and  communications.  These  missions 
will  require  stringent  pointing  accuracy,  shape  control,  and 
vibration  suppression.  To  satisfy  these  demanding  mission 
requirements,  the  LSS  will  almost  certainly  require  an 
active,  regulator-type  controller  with  multiple  sensors  and 
actuators  located  throughout  the  structure.1-4*  Furthermore, 
given  the  size  of  an  LSS,  there  will  be  a  large  set  of 
admissible  sensor  and  actuator  locations.  The  controls 
engineer  then  faces  the  problem  of  selecting  a  limited 
number  of  sensor  and  actuator  locations  to  "best  achieve" 
the  LSS  mission.  The  term  "best  achieve"  in  this  paper 
means  achieving  the  LSS  output  specifications  with  minimal 
actuator  power. 

The  solution  to  this  sensor  and  actuator  selection  (SAS) 
problem  needs  at  least  the  following  ingredients: 

(1)  A  specific  closed-loop  control  law  structure 

(2)  A  technique  for  systemically  adjusting  (tuning) 
control  law  parameters  to  achieve  output 
specifications  with  minimal  power 

(3)  A  technique  to  evaluate  the  effectiveness  of 
possible  sensor  and  actuator  configurations  in 
achieving  output  specifications  with  minimal  power 

This  paper  incorporates  the  above  ingredients  into  an 
algorithm  which  solves  the  SAS  problem  for  an  LSS  modeled  as 


a  linear  stochastic  system.  Some  necessary  background 
information  on  linear  stochastic  systems  and  linear 
quadratic  gaussian  (LQG)  control  theory  are  presented  in 
Section  2.0  along  with  a  formal  statement  of  the  SAS 
problem.  Section  3.0  contains  two  important  selection 
theorems,  four  pertinent  facts,  and  the  general  flow  of  the 
SAS  algorithm.  Results  from  the  SAS  algorithm  applied  to 
the  controller  design  for  a  large  space  telescope  are 
presented  in  Section  4.0,  and  Section  5.0  contains 
concluding  remarks. 


2 . 0  BACKGROUND 

The  SAS  algorithm  presented  in  this  paper  addresses 
controlling  an  LSS  described  by  the  following  linear 
stochastic  model: 

n  m  p 

x(t)  =  Ax(t)  +  Bu(t)  +  Dw(t)  ,  x€R  ,  u€R  ,  w€R 
x (t  )  =  x 

°  °  XXX 

D  *  [  B  D  ]  ,  w(t)  a  [  w  (t)  w  (t)  ] 

(2.1) 

y(t)  =  Cx(t),  y€R  (system  outputs) 

X 

z(t)  =  Mx(t)  +  v(t)  ,  z€R  (system  measurements) 
with  noise  characteristics 

X  0  0 

0  WS(t-T)  0 
0  0  VS(t-T) 


(2.2)  E 

X 

o 

w(t) 

r  X  ,  W  (t)  ,  V  (t)1  = 

v(t) 

L  °  J 

E  x  *  0,  E  w ( t )  =  0,  E  v  ( t ) 


0,  W>0,  V>0 
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The  notation  B  implies  a  real  vector  space  of  dimension  i. 
The  expectation  operator  is  denoted  by  E.  The  superscript  T 
represents  matrix  transposition.  The  Dirac  delta  function 
is  denoted  by  6,  and  W>0  implies  W  is  a  positive  definite 
matrix.  The  n-dimensional  vector  x(t)  represents  the  state 
of  the  system,  while  the  m-dimensional  vector  u(t)  contains 
actuator  signals.  The  regulated  system  outputs  are  defined 
by  the  k-dimens i onal  vector  y(t),  and  the  1-dimensional 
vector  z(t)  represents  the  system  measurements  (sensor 
information) .  The  white  noise  vector  process  w(t) 
represents  unmodeled  system  behavior  D  w  (t)  and  unmodeled 
actuator  behavior  B  w  (t) .  Unmodeled  sensor  behavior 
is  accounted  for  by  the  white  noise  vector  process  v(t). 

The  matrices  A,  B,  C,  D,  M,  W,  and  V  are  time-invariant  and 
appropriately  dimensioned.  The  matrix  B  has  no  zero 
columns,  and  the  matrices  C  and  M  have  no  zero  rows. 
Furthermore,  the  matrices  A,  B,  C,  D,  and  M  satisfy  the 
following  stabi 1 i zabi 1 ity  and  detectability  conditions: 

(2.3)  (A,B),  (A,D)  stabilizable ;  (A,C),  (A,M)  detectable. 

For  simplicity,  the  model  described  by  (2.1)-(2.3)  will  be 

identified  by  S(n,k,m,l).  The  arguments  n,  k,  m,  and  1 
represent  the  number  of  system  states,  outputs,  actuators, 
and  sensors  respectively. 

S(n,k,m,l)  specifications  usually  take  the  form 

2  2  2 

(2.4)  lim  Ey  (t)  =  E  y  L  a  ,  i  =  l,2,...k 

t  —  >—  1  a>  t  1 


(2.5)  lim  Eu  (t)  ■  E  u  i  n  ,  i  =  l,2,...m 

t->«  i  *  t  i 


i 

I 


i 


where  y  represents  the  ith  output  for  the  system,  and  u 
4  a  a  4 

represents  the  ith  input.  The  quantities  a  ,  |a  are 

i  ± 

mean  square  or  variance  constraints  on  the  ith  output  and 
input.  These  specifications  are  ideally  suited  for  the  LSS 
case  where  y  could  represent  a  line  of  sight  pointing 
error,  a  deflection  from  a  desired  shape,  etc.,  and  u 

X 

could  be  the  output  from  a  thruster,  torque  actuator,  etc. 

The  ingredients  of  the  SAS  solution  stated  in  the 
introduction  require  the  choice  of  a  specific  closed-loop 
control  law  which  is  readily  adjusted  to  meet  the  output 
specifications  with  minimal  power.  LQG  control  theory 
provides  such  a  law.  This  theory  originated  in  the  1960s 
with  the  foundational  work  of  Kalman,  Bucy ,  and  others. 
Since  that  time,  many  papers  and  texts  have  been  written 
further  clarifying  and  expanding  the  theory.  Refs.  7-9  are 
a  few  examples.  The  following  paragraphs  present  the 
fundamental  results  of  the  theory. 

For  a  system  S(n,k,m,l),  LQG  theory  guarantees  a  stable, 
steady-state,  closed-loop  control  system  minimizing 

T  T 

(2.6a)  V  =  E  (y  Q  y  +  u  R  u)  or 

M  a  n  3 

(2.6b)  V  =  2  E  y  q  +  2  Eur 

i  —1  »  1  1  1  _1  all 

where  q  and  r  are  the  ith  diagonal  entries  of  the 

i  t 

arbitrary,  diagonal,  and  positive  definite  weighting 
matrices  Q  and  R.  The  defining  equations  for  this  steady- 
state  controller  are 

—  IT  — 

(2.7a)  u (t )  =  -R  B  K  x(t)  - 


G  x  (t ) 


(2.7b) 


x(t)  =  A  x(t)  +  F  z(t)  ,  F  -  P  M  V 

c  ^ 

A  »  A  +  BG  -  FM,  x(t  )  =  x 


(2.7c)  KA  +  A  K  -  KBR  B  K  +  C  QC  =  0  (Steady  state 

Control  Riccati 
Equation) 

T  T  —  1  X 

(2.7d)  PA  +  AP  -  PM  V  MP  +  DWD  =  0  (Steady  state 

Filter  Riccati 
Equation) 


Substituting  this  controller  into  S(n,k,m,l),  the 
following  closed-loop  system  S(2n,  k+m,  m+1 ,  0)  results 

_  _  _  _  _  X  -T  X  _  XXX 

x=Ax+Bw,  x»(x  x  )  ,  w  »  (w  v) 


(2.8a) 


y=Cx,  V=yQy,  y«(y  u) 


(2 . 8b) 


A  BG 
FM  A 


.  B 


D  0 
0  F 


.  c 


0 


Q  0 
0  R 


Given  the  noise  characteristics  of  S(n,k,m,l),  the  steady 
state  variance  matrix  for  S(2n,  k+m,  m+1,  0)  is  known  to  be 

X+P  X 


(2.9)  E  <  (x  -  0)  (x  -  0)  >  -  E  <  x  x  >  = 


(2.10)  X  -  lim  X(t)  -  E  <  (x  -  0)  (x  -  0)  >  ■  E  {  x  x  > 

_  t  —  >—  v  m> 


where  X  is  the  solution  of  the  steady-state  Lyapunov 
equation 


X  —  X  —  1 

(2.11)  X(A  +  BG)  +  (A  +BG)  X  +  PM  V  MP  =  0 


Also,  coupling  the  definition  of  x  in  (2.8a)  with  (2.9) 
gives 
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iSi 


(2.12)  X  ■  lim  X(t)  aE{xx>  =  X  +  P 


The  following  facts,  documented  in  Refs.  10-14,  show  the 


LQG  control  structure  is  extremely  powerful  in  attacking  the 


SAS  problem: 


FI.  Analytical  expressions  for  E  y  and  E  u  exist  and 

■  1  a  i 


require  no  additional  major  calculations  beyond  those  of  the 


controller 


(2.13)  E  y  =  T  C  (P+X)  C  1  ,  i  =  l,...k,  and  n  implies 

•  *  L  J14  the  ith  diagonal 

2  f-  —  t  -I  element  of  the  matrix 

(2.14)  E  u  =  GXG  ,  i  =  l,  .  .  .m 

F2 .  A  successive  approximation  algorithm  (LQGWTS)  developed 


in  Refs.  10,  13,  and  14  can  be  used  to  adjust  the  elements 


of  Q  and  R  so  the  LQG  controller  achieves  the  output 


a  '  a 

specifications  E  y  =  cr  ,  i  ■  l,...k  with  minimal  power 

_  a  t  1 

m  St 

(i.e.,  2  E  u  =  min) .  LQGWTS  involves  no  gradient 

t-«  a  1 


calculations  or  complex  calculation  beyond  the  LQG 


controller  and  is  numerically  tractable  for  large  systems 


as  shown  in  Refs.  10,  13,  and  14. 


F3.  A  measure  developed  in  Refs.  10  and  11  identifies  the 


effectiveness  of  each  sensor  and  actuator  in  the  present  LQG 


controller  configuration.  The  mathematical  expressions  for 


these  effectiveness  values  are 


(2.15)  Y.  =  m  PLPm  V  (where  m  is  the  ith  column  of  M  ) 

i  iiii  i 


(2.16) 


•«=*  r  ~  t  i  r  t  a 

V  =  GXG  R  -  I  B  (K+L)  BW 


where  W  is  the  partition  of  W  corresponding  to  actuator 


noise  sources  (w  (t)),  and  L  satisfies  the  following: 


L  (A- PM  V  M)  +  (A- PM  V  M)  L  +  G  RG  =  0 


F2  and  F3  represent  recent  contributions  to  LQG  theory. 
Their  derivations  are  in  Ref.  10  but  are  omitted  here  for 
the  sake  of  brevity.  However,  Appendix  A  presents  a  brief 
summary  of  the  weight  selection  work,  and  Appendix  B 
provides  a  short  discussion  of  the  effectiveness  values. 

With  the  above  background,  the  SAS  problem  addressed  in 
this  paper  can  be  mathematically  stated: 

SAS  Problem  Statement 

Given :  A  system  (S,n,&,m,l)  which  has  only  m  out  of  m 

actuators  and  1  out  of  1  sensors  available  to 
design  a  steady  state  LQG  regulator  which 
must  achieve  a  set  of  output  variances  (cr*)  . 

Required :  Specify  the  closed-loop  system  which 

satisfies  the  following: 

m  a  22 

(2.17)  Min  2  E  u  subject  to  E  y  =  a  V  i=l , . . ,k 

1  —1  >  t  at  i 


3.0  SAS  ALGORITHM 


Before  presenting  the  algorithm,  two  selection  theorems 

and  some  pertinent  facts  should  be  stated. 

Theorem  1  Deletion  of  Noisy  Actuators 

For  a  system  S(n,k,m,l)  regulated  by  the  LQG 
controller  defined  in  (2.7),  the  deletion  of  an 
actuator  is  not  sufficient  for  V(m-l)>V(m)  where 

¥(m)  =  the  value  of  (2.6)  for  a  system  S(n,k,m,l) 
under  LQG  regulation 

Y(m-l)  *  the  value  of  (2.6)  for  a  system 

S(n,k,m-l,l)  'under  LQG  regulation 


The  proof  of  Theorem  1  is  given  in  Ref.  10.  More 

specifically,  Theorem  1  says  it  is  possible  for  an  LQG 

controller  to  do  better  (i.e.,  maintain  output 

specifications  at  a  lower  input  power)  when  a  noisy  actuator 

is  deleted.  This  means  the  desired  number  of  actuators  in  a 

control  system  might  be  less  than  the  allowable  (i.e.,  <  m) . 

Theorem  2  Deletion  of  Noisy  Sensors 

For  a  system  S(n,k,m,l)  under  the  regulation  of  the  LQG 
controller  defined  in  (2.7),  deletion  of  a  sensor 
cannot  reduce  the  V  of  (2 . 6)  . 

Theorem  2  clearly  says  that  noisy  sensors  never  degrade 

the  performance  of  an  LQG  controller.  The  proof  of  this 

theorem  is  contained  in  Refs.  10,  12,  and  15.  However, 

intuitively,  it  makes  sense  because  the  LQG  controller  uses 

a  Kalman  filter,  and  if  a  sensor  is  too  noisy  the  Kalman 

filter  essentially  ignores  the  information  from  that  sensor. 

Therefore,  Theorem  2  indicates  the  number  of  sensors  in  the 

controller  design  should  be  1  (i.e.,  the  allowable  number). 

The  following  facts  are  important  to  the  algorithm: 

act 

F4 .  Data  in  Refs.  10  and  16  show  relative  V  and  V 

&  1 

rankings  can  change  as  functions  of  Q  and  R. 

F5 .  Deleting  sensors  and  actuators  changes  the  Q  and  R 
needed  to  achieve  output  specifications  at  minimal  power. 

F6 .  The  results  from  Refs.  10,  11,  15,  and  17  show  the 
SAS  problem  is  highly  coupled  and  simultaneous  SAf%  works 
better  than  sequential  selection. 


F7  .  V  and  V  do  not  predict  the  loss  of  measurability  or 

X  * 

controllability . 

With  these  theorems  and  facts  stated,  the  general  flow  of 
the  SAS  algorithm  can  be  outlined. 


SAS  Algorithm  (General  Flow) 

1.  Given  the  system  S(n,k,m,l)  specify  the  number  of 
sensors  1  (1<1)  and  actuators  m  (m<m)  that  can  be 
used  in  the  controller  design. 

2.  Run  the  subroutine  LQGWTS  for  the  full  system 
S(n,k,m,l).  (LQGWTS  selects  the  weighting  matrices 
Q  and  R  so  the  LQG  controller  enables  the  closed- 
loop  system  to  satisfy  (2.17).) 

3.  With  the  resulting  closed-loop  system,  calculate 
V~"  and  V*44*  for  each  sensor  and  actuator  using 
(2.15)  and  (2.16)  . 

4.  Rank  the  actuators  and  sensors  according  to  the 
algebraic  value  of  V and  V**5*  . 

5.  Modify  the  current  system  to  eliminate  the  lowest 
ranking  sensor  and  actuator  and  all  others  of 
"nearly"  the  same  ranking,  if  measurability  and/or 
controllability  of  the  system  are  not  disturbed. 

6.  If  the  appropriate  number  of  actuators  and  sensors 
are  deleted  (i.e.,  m-m  and  1-1)  stop.  However,  in 
light  of  Theorem  1,  fewer  actuators  than  m  may  be 
better,  and  actuator  deletion  should  continue  until 

*•  a 

(3.1)  (1/m-)  2  E  u 

x  — i  —  * 

is  no  longer  decreasing.  (Note,  m*  is  the  number 
of  actuators  currently  in  the  design  (m-  <  m) ,  and 
Equation  3.1  is  the  average  power  per  actuator.) 

7.  Run  LQGWTS  on  the  reduced  system  and  return  to  3. 


F4  and  F5  dictate  the  order  of  Steps  2,  3,  and  4  in  the 
above  general  flow.  More  specifically,  since  the  ranking  of 
effectiveness  values  can  vary  as  a  function  of  Q  and  R,  the 


effectiveness  values  should  be  based  on  the  Q  and  R  that 


produce  a  closed-loop  system  satisfying  (2.17). 

Step  5  of  the  general  flow  results  from  facts  F6  and  F7 . 
Caution  should  be  exercised  in  deleting  sensors  and 
actuators  with  "nearly"  the  same  effectiveness  values  since 
F4  and  F5  indicate  any  perturbation  in  the  sensor  and 
actuator  configuration  can  alter  relative  effectiveness 
rankings  of  the  sensors  and  actuators.  Also,  the 
requirement  for  checking  measurability  and  controllability 
is  severe.  However,  as  a  result  of  the  work  by  Skelton  and 
Hughes  in  Ref.  18,  the  controllability  and  measurability  of 
systems  in  modal  form  (i.e.,  the  form  of  the  telescope  model 
in  Appendix  C)  can  be  done  by  inspection  for  systems  of  any 
dimension. 

Theorems  1  and  2  provide  the  rationale  for  Step  6.  If  1 
sensors  are  allowed,  then  according  to  Theorem  2,  1  sensors 
should  be  used.  Theorem  1  invalidates  this  philosophy  for 
actuators . 

4.0  SPACE  TELESCOPE  EXAMPLE 

This  section  presents  the  results  of  the  SAS  algorithm 
when  applied  to  the  model  of  a  solar  optical  telescope. 
Before  discussing  the  results,  a  brief  model  description  is 
presented . 

The  solar  optical  telescope  model  was  developed  in  1980  by 
the  Charles  Stark  Draper  Laboratory  (CSDL)  primarily  for 
providing  a  minimum  complexity  structure  to  evaluate  LSS 


control  design  techniques.  The  outputs  for  the  model  are 
the  telescope  line  of  sight  angles  about  the  x  and  y  axes 
(see  Figure  4.0)  and  the  focal  length  (defocus)  of  the 
lenses  located  at  the  top  and  bottom  of  the  telescope. 


The  12  nodes  shown  in  Figure  4.0  represent  the  spatial 
locations  for  45  admissible  sensors  and  21  admissible 
actuators.  Tables  4.0  -  4.2  provide  specific  type, 
location,  direction,  and  label  information  for  the 
admissible  sensors  and  actuators  and  also  a  listing  of  th 
output  specifications  (a) . 
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Table  4.0  Telescope 

sensor 

description 

Sensor 

#  Tvpe _ Nodal  Location _ Direction 

Label 

1 

Line  of  Sight  Angle 

_ _ 

X 

LOS* 

2 

1  1 

— 

Y 

LOSv 

3 

Defocus 

— 

— 

Defocus 

4 

Linear  Displacement 

1 

Y 

Y1 

5 

1 1 

1 

Z 

Z1 

6 

1  I 

2 

Z 

Z2 

7 

1  1 

3 

X 

X3 

8 

1  1 

3 

Y 

Y3 

9 

1  1 

3 

Z 

Z3 

10 

1  1 

4 

Z 

Z4 

11 

1  1 

5 

X 

X5 

12 

1  1 

5 

Y 

Y5 

13 

1  1 

5 

Z 

Z5 

14 

I  1 

6 

Z 

Z6 

15 

1  1 

7 

Y 

Y7 

16 

1 1 

7 

Z 

17 

17 

1  > 

8 

Z 

18 

18 

1  1 

9 

Z 

19 

19 

10 

Z 

Z10 

20 

1  1 

11 

X 

Xll 

21 

11 

Y 

Yll 

22 

11 

Z 

Zll 

23 

1  1 

12 

Y 

Y12 

24 

II 

12 

Z 

Z12 

25 

Linear  Rate 

1 

Y 

LRY1 

26 

1  1 

1 

Z 

LRZ1 

27 

1  1 

2 

Z 

LRZ2 

28 

1  1 

3 

X 

LRX3 

29 

1 1 

3 

Y 

LRY3 

30 

1  1 

3 

Z 

LRZ3 

31 

l< 

4 

Z 

LRZ4 

32 

1  1 

5 

X 

LRX5 

33 

1  1 

5 

Y 

LRY5 

34 

1  1 

5 

Z 

LRZ5 

35 

1  1 

6 

Z 

LRZ6 

36 

3  1 

7 

Y 

LRY7 

37 

1  1 

7 

Z 

LRZ7 

38 

1  1 

8 

Z 

LRZ8 

39 

1  1 

9 

Z 

LRZ9 

40 

1  1 

10 

Z 

LRZ10 

41 

1  1 

11 

X 

LRX11 

42 

1  1 

11 

Y 

LRY11 

43 

1 1 

11 

Z 

LRZ11 

44 

1  1 

12 

Y 

LRY12 

45 

1  1 

12 

Z 

LRZ12 

Actuator  (Force 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Nodal  Location 


1 

1 

2 

3 

3 

3 

4 

5 
5 

5 

6 
7 

7 

8 

9 

10 
11 
11 
11 
12 
12 


Direction 


Y 
Z 
Z 
X 

Y 
Z 
Z 
X 

Y 
Z 
Z 

Y 
Z 
Z 
Z 
Z 
X 

Y 
Z 

Y 
Z 


Label 


Table  4.2  Telescope  specifications 


Optical  line  of  sight 

I  I 

Defocus 


Label _ Value. 

LOS*,  65.2  sec 

LOS^ 

Defocus  0.001  mm 


CSDL  used  the  program  NASTRAN  to  develop  data  for  the 
first  44  mode  shapes  of  the  telescope  and  provide  location 
information  for  the  two  sinusoidal  disturbances  S*  and  Sa 
shown  in  Figure  4.0.  The  final  design  model  chosen  for  this 
example  was  a  10  mode,  20  state,  linear  stochastic  model 
coupled  with  a  2  mode,  4  state  linear  stochastic  model  of 


the  disturbances  (i.e.,  a  model  St«.i  •  (24 , 3 ,21 , 45)  )  .  The 


wi  wi’w.uuur^jumjui  ^ mmm^hm w  tfi yn ininnn vn  wv¥  ■.■vi^i.’vv^^  iva'.v  v*^r 
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development  of  this  model  from  the  NASTRAN  data  is  presented 
in  Ref.  10.  The  model  data  is  in  Appendix  C. 

The  SAS  problem  for  the  telescope  is 
Telescope  SAS  Problem 

Given :  Sc.i . <24 , 3 , 2 1 , 45)  with  only  12  actuators  and  12 

sensors  available  for  designing  an  LQG  regulator  to 
achieve  the  ( a )  of  Table  4.2. 

Required :  Specify  the  closed-loop  system  which 

satisfies  (2.17). 

Table  4.3  shows  the  iterations  of  the  SAS  algorithm  when 
applied  to  the  space  telescope.  Three  observations  are 
worth  noting.  First,  in  Iteration  #1  of  the  algorithm,  the 
measurements  of  the  outputs  were  deleted.  This  might  appear 
disturbing,  but  two  things  should  be  remembered:  the 
algorithm  guarantees  the  system  remains  measureable  and  the 
goal  of  the  Kalman-Bucy  filter  is  to  estimate  the  state  of 
the  system  and  not  necessarily  the  output.  In  the  telescope 
example,  the  noise  present  on  the  output  measurements  made 
them  unattractive  to  the  filter.  Second,  despite  sensor 
deletion,  the  average  input  value  decreased  from  Iterations 
1-3.  This  is  dramatic  empirical  proof  of  Theorem  1  and 
indicates  the  controller  is  more  efficient  at  holding  the 
system  outputs  to  specification  with  19  sensors  than  it  is 
with  20  or  21.  Furthermore,  the  rise  of  average  actuator 
power  at  Iteration  4  suggests  the  optimal  number  of 
actuators  for  this  design  lies  between  19  and  15.  Third, 
the  algorithm  converged  to  a  configuration  in  only  10 
iterations.  A  direct  search  algorithm  could  have  taken 


thousands  of  iterations  even  with  the  aid  of  gradient 
calculat ions. 


Table  4.3  Telescope  SAS  algorithm  results 


Y3(.  008197) 
Y1 (.008195) 
L0SY(.  000015) 
LOSXC. 000002) 
DEFOCUS  (0) 


FX11  (-7.980) 


7.483 


45/21 


Z5  (  .02955) 
Z7(. 02955) 


FX5  (-8.796) 


7 . 393 


40/20 


Y12  ( 
Y7( 
Yll  ( 
YS  ( 
Z2  ( 
Z4( 
Zl( 
Z3( 
Z6( 
Z8( 


.03848) 

.03848) 

.03848) 

.03484) 

.03835) 

.03835) 

.03570) 

.03570) 

.03460) 

.08460) 


FY7(-4 .003) 
FY12  (-4.003) 
FY11  (-4.003) 
FY5  (-4.003) 


7.221 


38/19 


Z10(. 08602) 
Z9(. 08602) 
X3(. 08277) 


FZ11 (-4 .8877) 
FZ12 (-4.8877) 


7.789 


28/15 


LRY3(. 10815) 
LRY1 (.10813) 


FZ9  (-5 .414) 


8.155 


25/13 


LRZ7(.  12273) 
LRZ5  (  .11553) 


8 . 339 


23/12 


LRZ1 (.16805) 
Zll  (.16129) 
LRZ6(. 13838) 
LRZ8(. 13615) 


8.359 


21/12 


LRX3 ( . 2237) 
LRZ2  (.2191) 
LRZ4  (.2077) 
LRZ3  (  . 1960) 


8.416 


17/12 


LRZ9(. 47256) 


8 .492 


13/12 


8.535 


12/12 


configurations  have  been  tested  and  none  performed  as  well 
as  this  configuration.10 

Table  4.4  shows  the  output  specifications  achieved  by  this 
regulator  and  the  root-mean-square  (rms)  power  of  each 
actuat  or . 


Table  4.4  Telescope  performance 


Output  # 

<  E— y  i 3 

Actuator  # 

<  E-ui  2 

(minimum  achievable) 

1 (LOSx) 

65.2  sec 

1 (FY1) 

. 0244N 

2 (LOSy) 

1  1 

2 (FZ1) 

. 0303N 

3 (DEFOCUS) 

. 0002mm 

3 (FZ2) 

. 0297N 

4  (FX3) 

. 0324N 

5  (FY3) 

. 025 IN 

6  (FZ3) 

. 0265N 

7 (FZ4) 

. 0338N 

8 (FZ5) 

. 0245N 

9 (FZ6) 

. 0242N 

10 (FZ7) 

.  0221N 

11  (FZ8) 

. 0273N 

12  (FZ10) 

. 0437N 

The  rms  power,  expressed  in  Newtons  (N) ,  gives  the 
controls  engineer  actuator  sizing  information  for  the 
controller  design.  It  also  provides  a  measure  of  the 
importance  of  each  actuator  in  the  controller.  This 
information  is  useful  in  reliability  and  redundancy 
considerations  for  the  actuator  configuration. 

A  final  note  is  the  value  of  the  defocus  output.  The 
value  (.0002mm)  is  well  below  the  u  of  ,001mm.  This  seems 
to  indicate  the  controller  is  working  too  hard  and  less 
actuator  power  would  be  required  if  the  defocus  operated  at 
.001mm.  However,  for  this  model,  the  defocus  output  is 


■  A.— V.w.^ 


linearly  dependent  upon  LOS*,  and  LOSv  (i.e.,  the  defocus 
output  could  have  been  deleted  from  the  model) .  The 
algorithm  LQGWTS  looks  for  output  dependency  and  assigns 
dependant  outputs  a  zero  weighting  to  delete  them  from  the 
model . 

5.0  CONCLUSION 

In  summary,  the  SAS  algorithm  aids  the  controls  engineer 
in  specifying  a  sensor  and  actuator  configuration  for 
regulation  of  large  scale,  linear,  stochastic  systems.  The 
algorithm  uses  an  LQG  controller,  an  efficient  weight 
selection  technique  based  upon  successive  approximation,  and 
a  measure  of  sensor  and  actuator  effectiveness  to  specify  a 
final  sensor  and  actuator  configuration.  This  configuration 
enables  the  closed-loop  system  to  meet  output  specifications 
with  minimal  input  power.  Also,  the  algorithm  involves  no 
complex  gradient  calculations  and  has  proven  numerically 
tractable  for  large  linear  models .  This  tractability  is 
demonstrated  by  the  space  telescope  example  presented  in 
Section  4.0  and  the  hoop-column  antenna  example  of  Refs.  10 
and  12.  Finally,  the  algorithm  provides  the  controls 
engineer  information  on  the  important  design  issues  of 
actuator  sizing,  reliability,  redundancy,  and  optimal 


number . 
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APPENDIX  A  LQG  WEIGHT  SELECTION 


The  successive  approximation  technique  (LQGWTS)  adjusts 
the  diagonal  Q  and  R  matrices  in  the  LQG  controller  so  the 
closed-loop  system  satisfies  variance  constraints  on  inputs 
outputs,  or  both.  The  technique  evolved  through  Refs.  14, 
13,  and  10.  This  appendix  presents  a  summary  of  the  LQGWTS 
theory  applying  to  the  minimum  power  SAS  problem. 

For  the  minimum  power  problem,  each  iteration  of  the  SAS 
algorithm  requires  LQGWTS  to  design  a  closed-loop  system 

meeting  output  variance  specifications,  while  minimizing 

m  a 

actuator  power  (2  E  u  )  .  LQGWTS  does  this  by  first  setting 

*.  —  i 

R  equal  to  the  identity  matrix.  This  allows  the  LQG 

controller  to  consider  all  actuators  equally  and  permits  th 

m  a 

term  2  E  u  to  appear  directly  in  the  cost  function  V  of 

*  .t  —  i 

(2.6),  thus  assuring  its  minimization.  LQGWTS  then 
adjusts  Q  using  the  following  iterative  equations: 

q  (1)  =  1/a  ,  i  =  1 , . . ,k 

±  & 

(A .  1)  q  <j  +  l)  =  [  E  y  <j)/a  ]  •  q  (j)  ,  j  =  1,2,... 

i  «ii  i 

pwr(l)  =  1,  pwr(j)  =  j-1  V  j>l 
where  q  is  the  ith  diagonal  entry  of  Q,  and  j  is  the  jth 
iteration  of  LQGWTS. 

The  update  expression  in  (A.l)  always  adjusts  q  in  the 
right  direction.  The  exponent  pwr(j)  is  used  to  speed  up 
convergence  (i.e.,  increase  step  size).  To  prevent 


overshoot  LQGWTS  incorporates  a  descent  function  defined  as 

a  a 

(A. 2)  descent  (j  +  1)  =  max  [  (E  y  (j  +  l)/o  -1),  0] 

si  i 

If  descent  (j  +  1)  >  descent  (j),  LQGWTS  returns  to  the 

jth  iteration,  resets  the  pwr  function  to  1  (i.e.,  the 
smallest  step  size),  and  continues. 

As  mentioned  in  Section  4.0,  LQGWTS  checks  for  dependent 
outputs  and  removes  them  from  the  model .  The  check  LQGWTS 
applies  is 

(A. 3)  If  <  €/cti  then  qi»0  (€<<1) 

where  setting  qi=0  removes  the  ith  output  from  the  model. 

When  used  for  the  telescope  SAS  example,  LQGWTS  converged 
in  9  or  less  iterations  for  each  iteration  of  the  SAS 
algorithm.  The  major  advantage  of  LQGWTS  is  its  lack  of 
cumbersome  gradient  calculations.  Additionally,  it  has 
demonstrated  "good"  convergence  properties  on  substantial 
problems  in  Refs.  10,  12,  13,  and  14. 


APPENDIX  B  ACTUATOR  AND  SENSOR  EFFECTIVENESS  VALUES 


The  effectiveness  values  in  (2.15)  and  (2.16)  are  derived 
from  closed-loop  input  cost  analysis  (CICA)  and  closed-loop 
output  cost  analysis  (COCA) .  Both  CICA  and  COCA  are 
extensions  of  component  cost  analysis  (CCA)  developed  by 
Skelton  and  co-workers  in  the  late  seventies  and  early 
eighties.  CICA  and  COCA,  along  with  the  sensor  and  actuator 
effectiveness  values,  are  discussed  in  detail  in  Refs.  10- 
12.  This  appendix  briefly  presents  the  important  details 
applying  to  the  telescope  SAS  problem. 

The  fundamental  concept  behind  CICA  and  COCA  is  to 
identify  the  contribution  each  input  and  output  of  (2.8)  is 
making  to  the  cost  functional  V  of  (2.6).  With  these 
contributions  identified,  they  are  combined  to  determine  the 
effectiveness  of  each  input  and  output  to  the  control 
effort . 

Using  the  closed-loop  model  (2.8),  the  following  input  and 
output  costs  are  defined: 

u  _t  _  _ 

(B.l)  V  *  1/2E  {  (a?y  Q  y/c?u  )u  >  (i.e.,  contributi  on  ith 

*  “  11  control  is  making  to  V) 

w  _ T  _  _ 

(B.2)  V  “  ■  1/2E  <  (o?y  Q  y/c?w«  )w«  >  (i.e.,  contribution  ith 
1  “  14  actuator  noise  is 

making  to  V) 

V  _T_  _ 

(B.3)  V  *  1/2E  <  (oly  Q  y/ci^v  )v  >  (i.e.,  contribut  ion  ith 

1  *  11  sensor  noise  is  making 

to  V) 

Applying  these  definitions  to  the  closed-loop  model  (2.8) 


yields  the  following  formulas: 


CB.4)  V 


[GXG  R]  ,  i  =  1 ,  .  .  ,  .m 

41  (u  subscript  means  ith  diagonal 
element  of  matrix) 


(B  .5)  V  •  ■ 


[B  CK+L) BW  ]  ,  i  =  1, 

1 1 


,m 


(B . 6)  X  *  [F  LFV]  ,  i  =  1, _ 1 

&  ii 

where  the  matrix  L  satisfies  the  following  steady  state 
Lyapunov  equation: 

X  X 

(B.7)  L(A-FM)  +  (A-FM)  L  +  G  RG  =  0 

Equations  (B.4)-(B.6)  and  some  useful  properties  are 
derived  in  Refs.  10  and  11.  These  equations  are  the  basis 

act  ««n 

for  V  and  V  defined  in  (2.15)  and  (2.16). 

&  ft 

u 

As  noted  above,  V  represents  the  contribution  that  u  is 


making  to  V.  Since  the  function  of  LQG  theory  is  to  use  u 


to  minimize  V  (for  a  given  Q  and  R)  ,  a  large  \L  means  u  is 

i  i 


important  to  the  minimization  effort.  Alternatively,  the 

M 

contribution  the  ith  actuator  noise  source  makes  to  V  (V  *) 

ft 

is  clearly  an  undesirable  contribution.  These  facts  coupled 


with  additional  rationale  discussed  in  Refs.  10-13  produce 


ACt  U  M 

(B .  8)  V  -  V  -  V  - 

ft  ft  ft 


Substituting  (B.4)  and  (B.5)  into  (B.8)  gives  (2.16) . 

A  modification  to  (B.8)  accounting  for  constrained  actuator 
power  is  presented  in  Ref.  10.  A  modification  is  also 
developed  in  Ref.  19  to  account  for  noise  correlations. 

As  previously  noted,  V  is  the  contribution  the  ith 

i 

sensor  noise  source  makes  to  V.  Therefore,  those  sensors 


with  larger  values  for  V  are  contributing  more  noise  to  V 


than  those  with  a  smaller  V  and,  at  first  glance,  appear  to 
be  candidates  for  deletion.  However,  the  sensor 
measurements  for  systems  S(n,k,m,l)  under  LQG  regulation  are 
passed  through  a  Kalman-Bucy  filter.  One  of  the  purposes  of 
this  filter  is  to  de-emphasize  or  throw  out  measurements 
which  have  more  noise  than  estimation  information. 

Therefore,  any  noise  source  making  a  "large"  contribution  to 
V  comes  from  a  sensor  which  is  making  an  even  larger 
contribution  to  the  estimation  information  necessary  to 
minimize  V,  Consequently,  the  following  effectiveness  value 
is  defined: 


CB.9)  V 


wry  x 


Substituting  (B.6)  and  the  definition  for  F  in  (2.8)  into 
(B.9)  yields  (2.15).  Additionally,  Ref.  19  presents  a 
modification  to  (B.9)  accounting  for  noise  correlations. 


APPENDIX  C  TELESCOPE  MODEL  DATA 


This  appendix  presents  the  contents  of  the  A,  B,  C,  D,  M, 
W,  and  V  matrices  for  the  telescope  model  S  (24 , 3 , 21 , 45)  . 


(C.l)  A  = 


£2  S  0to 

0. 

0m  Oim 

0.  S2.  0xo 

,  B  = 

BE 

,  D  = 

BE  la 

CU  0  4> 

0. 

0.  0a 

- 

BR 

BR  0a 

P 

(C .  2)  C  =  [  CE 


CR  02]  ,  M  = 


CE  O12  CR  CP 
ME  O12  MR  CP 
CP  ME  CP  MR 


(C .  3) 


r  w- 

Oai 

Vx 

0a 

0a 

L  O' 

W° 

,  v  = 

Oaa 
-  Oa  l 

Va 

Oat 

Oaa 

Va 

The  terms  in  the  above  matrices  are  defined  below: 


Q 


[- 


Om  I. 

w2  -2 a 


0.  *  null  matrix  with  8  rows 
I«  ■  8x8  identity  matrix 
w2  «  8x8  diagonal  matrix  of  squared 
modal  frequencies 

2cx  *  8x8  diagonal  modal  damping  matrix 


=  diag [ . 835 ,  2.736,  3.971,  4.378,  7.746,  13.175, 
13.339,  59.112]  (rad2/sec2) 


2a  =  diag[. 00183,  .00331,  .00399,  .00418,  .00557, 

.00726,  .00730,  .0154]  (rad/sec) 


8x2  coupling  matrix  (defined  in  table  C . 0) 
w.2 


0a  la 

2-  =  -w.2  -2a. 


■  2x2  diagonal  matrix  of 
squared  disturbance 
frequencies 

2a.  ■  2x2  diagonal  disturbance 
damping  matrix 


w«2  =  diag[3947.8,  986.96]  (rad2/sec2) 
2a.  *  diag[.1257,  .0628]  (rad/sec) 


« 


0a  la 
0a  0a 


* 


telescope  rigid  body  rotational 
modes  (9m ,0v) 


8x21  matrix  transpose  of  modal  eigenvectors 
(defined  in  table  C.l) 

2x21  matrix  transpose  of  rigid  body 
eigenvectors  (defined  in  table  C.2) 

3x8  matrix  of  linear  combinations  of  modal 
eigenvectors  required  by  system  output  (defined 
in  table  C  .  3) 

«  null  matrix  with  12  columns 

3x2  matrix  of  rigid  body  contribution  to  system 
outputs  (defined  in  table  C.4) 

21x8  matrix  of  linear  combinations  of  modal 
eigenvectors  required  by  system  measurements 
(defined  in  table  C.5) 

21x2  matrix  of  rigid  body  contribution  to  system 
measurements  (defined  in  table  C . 6) 

(0.1) lax  (N2)  ,  W  ■  (3 . 95)  la  (N2)  , 

(.0001) Ia  (rad2),  V«  «  (.000001) Iaa  (m2) , 

(1x10— 7) lax  (m2/s2) 


Table  C.O  The  S  matrix  (8x2) 


-1 .5105E-03 
8 . 0569E-04 
-5 . 1719E-03 
1 .0263E-03 
7 . 8788E-03 
3 . 0956E-03 
3 . 5394E-03 
1 .3823E-04 


6 . 1789E-03 
-9 . 4151E-05 
-7 . 4586E-04 
1 .4033E-03 
-5 . 0280E-03 
-4 .2834E-03 
-5 . 3071E-04 
2 . 214  IE-02 


Table  C.l  The  BE  matrix  (8x21) 


ol:  1 

2 

3 

4 

5 

6 . 740E-05 

-1 .520E-03 

-1 .507E-03 

1 .946E-03 

-6 . 628E-05 

4 . 615E-03 

9 . 353E-04 

7 . 917E-04 

-7 . 844E-03 

4 . 612E-03 

2 . 582E-04 

-5 . 140E-03 

-5 . 175E-03 

2 . 359E-06 

-2 . 581E-04 

4 . 927E-03 

-7 . 887E-04 

1 .479E-03 

1 .084E-06 

-4 . 930E-03 

2 . 786E-04 

7 . 798E-03 

7 . 885E-03 

-7 . 174E-06 

2 .  784E-04 

1 .557E-02 

2.626E-03 

-4 .516E-03 

-6 . 149E-05 

1 .558E-02 

3 . 792E-04 

-3 . 451E-03 

-3. 544 E- 03 

1 .705E-02 

4 . 814E-04 

6 . 065E-05 

8.11 3E-05 

1 .475E-04 

4 . 850E-04 

6 . 065E-05 

w\sjo1j _ 6_ 

7 

8 

9 

10 

1 .521E-03 

1 .507E-03 

-6 . 035E-03 

8 . 661E-05 

-1 . 5 17E-03 

-9 . 354E-04 

-7 . 91 IE-04 

-1 .575E-03 

-4 . 944E-03 

9 . 175E-04 

-5 . 140E-03 

-5 . 175E-03 

-3 . 486E-08 

-3 . 288E-04 

-5 . 152E-03 

-7 . 889E-04 

1 .480E-03 

-4 . 4 65E-09 

-3.893E-04 

-3 . 367E-04 

7 . 798E-03 

7 . 885E-03 

1 .040E-07 

4 . 498E-04 

7 . 831E-03 

2 . 607E-03 

-4 . 54 6E-03 

2 . 331E-06 

1 .257E-03 

1 . 212E-03 

3 . 468E-03 

3 . 5 14E-03 

-4 . 789E-04 

-5 . 916E-04 

-3 . 473E-03 

-8 . 125E-05 

-1 .471E-04 

-3 . 086E-04 

6 . 402E-05 

9 . 512E-05 

w\col :  11 

12 

13 

14 

15 

-1 .510E-03 

-8 . 696E-05 

1 .517E-03 

1 .510E-03 

-1 .522E-03 

8.049E-04 

4 . 942E-03 

-9 . 180E-04 

-8 . 050E-04 

1 .042E-03 

-5 . 172E-03 

-3 . 290E-04 

-5 . 152E-03 

-5 . 172E-03 

-5 . 132E-03 

1 .026E-03 

-3 . 896E-04 

-3 . 366E-04 

1 .027E-03 

-1 .928E-03 

7 . 879E-03 

4 . 502E-04 

7 . 831E-03 

7 . 879E-03 

7 . 786E-03 

-3 . 095E-03 

1 .250E-03 

1 . 190E-03 

-3 . 124E-03 

6 . 250E-03 

-3 . 540E-03 

5 . 986E-04 

3 . 4  81E-03 

3 . 519E-03 

-3 . 405E-03 

1 .417E-04 

-6 . 462E-05 

-9 . 552E-05 

-1 .408E-04 

6 . 075E-05 

w\col !  16 

17 

18 

19 

20 

1 .522E-03 

-6.313E-03 

8 . 542E-05 

-1 .504E-03 

-8 . 536E-05 

-1 .043E-03 

1 .452E-02 

-4 . 945E-03 

6 . 931E-04 

4 . 945E-03 

-5 . 132E-03 

2 . 594E-07 

-3 . 290E-04 

-5.195E-03 

-3.290E-04 

-1 .929E-03 

-3. 986E-07 

-3 . 892E-04 

2 . 617E-03 

-3 . 894E-04 

7 . 786E-03 

-5 . 702E-07 

4 . 5 01E-04 

7 . 933E-03 

4 . 504E-04 

6 . 236E-03 

-1 .038E-05 

1 .256E-03 

-8. 125 E- 03 

1  . 24  8E-03 

3 . 446E-03 

1 .473E-03 

-5 , 928E-04 

-3. 611 E- 03 

6 . 008E-04 

-6 . 123E-05 

-5 . 088E-04 

5 . 993E-05 

1 .824E-04 

-5 . 937E-05 

w\col :  21 

1 .503E-03 

-6 . 904E-04 

-5 . 195E-03 

2 . 618E-03 

7 . 933E-03 

-8 . 162E-03 

3.558E-03 

-1 .824E-04 

S1 


Table  C.5  The  ME  matrix  (21x8) 


6 . 740E-Q5 
1 .520E-03 
1 .507E-03 
1 .946E-03 
6.628E-05 
1 .521E-03 
1 . 507E-03 
6 . 035E-03 
8 . 661E-05 
1 . 517E-03 
1 .510E-03 
8 . 696E-05 
1 . 5 17E-03 
1 .510E-03 
1 .522E-03 
1 .522E-03 
6.313E-03 
8 . 542E-05 
1 .504E-03 
8 . 536E-05 
1 .503E-03 


■4 . 615E-03 
9 . 353E-04 
7 . 917E-04 
-7 . 844E-03 
4 . 612E-03 
-9 . 354E-04 
-7.911E-04 
-1 .575E-03 
■4 . 944E-03 
9 . 175E-04 
8 . 04  9E-04 
4 . 942E-03 
-9 . 180E-04 
-8 . 050E-04 
1 . 042E-03 
-1 . 04  3E-03 
1 .452E-02 
-4 . 945E-03 
6 . 931E-04 
4 . 945E-03 
-6 . 904E-04 


-2 . 5  82E-04 
-5 . 140E-03 
-5 . 175E-03 
2 . 359E-06 
-2 . 581E-04 
-5 . 140E-03 
-5 . 175E-03 
-3 . 4  86E-08 
-3 . 2  88E-04 
-5 . 152E-03 
-5 . 172E-03 
-3 . 290E-04 
-5 . 152E-03 
-5 . 172E-03 
-5 . 1 32E-03 
-5 . 1 32E-03 
2 . 5  94E-07 
-3 . 290E-04 
-5 . 195E-03 
-3 . 290E-04 
-5 . 1 95E-03 


-4 . 927E-03 
-7 . 887E-04 
1 .479E-03 
1 .083E-06 
-4 . 930E-03 
-7 . 889E-04 
1 .480E-03 
-4 . 4  65E-09 
-3 . 893E-04 
-3 . 367E-04 
1 . 026E-03 
-3 . 896E-04 
-3 . 366E-04 
1 .027E-03 
-1 .928E-03 
-1 .929E-03 
-3 . 986E-07 
-3 . 892E-04 
2 . 617E-03 
-3 . 894E-04 
2 . 61 8E-03 


1 .557E-02 
2 . 626E-03 
-4 . 5 16E-03 
-6 . 149E-05 
1 .558E-02 
2 . 607E-03 
-4 . 546E-03 
2 . 331E-06 
1 .257E-03 
0  1 . 212E-03 

1  -3 . 095E-03 

2  1 . 250E-03 

3  1.190E-03 
3 . 124E-03 
6 . 250E-03 
6 . 2  36E-03 
1 .037E-05 
1 .256E-03 
8 . 125E-03 
1 .248E-03 
8 . 162E-03 


-3 . 792E-04 
-3 . 451E-03 
-3 . 544E-03 
1 .705E-02 
4 . 814E-04 
3 . 468E-03 
3 . 5 15E-03 
-4 , 789E-04 
-5 . 916E-04 
-3 , 473E-03 
-3 . 540E-03 
5 . 986E-04 
3 , 481E-03 
3 .519E-03 
-3 . 405E-03 
3.446E-03 
1 .473E-03 
-5 . 928E-04 
-3.611E-03 
6 . 008E-04 
3 . 558E-03 


-6 . 065E-05 
8.11 3E-05 
1 .475E-04 
4 . 850E-04 
6 . 065E-05 
-8 . 125E-05 
-1 .471E-04 
-3 . 086E-04 
6 . 4  02E-05 
9.512E-05 
1 . 4 17E-04 
-6 . 462E-05 
-9 . 552E-05 
-1 .408E-04 
6 . 075E-05 
-6. 123E-05 
-5 . 088E-04 
5 . 993E-05 
1 .824E-04 
-5 . 937E-05 
-1 .824E-04 


2 . 78~E-04 
7 . 798E-03 
7 . 885E-03 
-7 . 174E-06 
2 . 784E-04 
7 . 798E-03 
7 . 885E-03 
1 . 040E-07 
4 . 4  98E-04 
7 . 831E-03 
7 . 879E-03 
4 . 502E-04 
7 . 831E-03 
7 . 879E-03 
7 . 786E-03 
7 . 786E-03 
-5 . 702E-07 
4 . 501E-04 
7 . 933E-03 
4 . 504E-04 
7 . 933E-03 
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Table  C.6  The  MR  matrix  (21x2) 


row\col : _ 1 

2 

1 

5 . 632E+00 

-4 . 545E-07 

2 

4 . 463E+00 

4 . 000E+QQ 

3 

-5 . 537E+00 

4 . OOOE+OO 

4 

0 

-5 . 632E+00 

5 

5 , 632E+00 

4 . 545E-07 

6 

4  .463E+00 

-4 .000E+00 

7 

-5 .  537E+00 

-4 .000E+00 

8 

0 

1 .437E+01 

9 

-1 .437E+01 

-4 . 545E-07 

10 

2 . 463E+00 

4 .000E+00 

11 

-3 . 537E+00 

4 .000E+00 

12 

-1 .437E+01 

4 . 545E-07 

13 

2  .463E+00 

-4 .OOOE+OO 

14 

-3 . 537E+00 

-4 .OOOE+OO 

15 

9 . 463E+00 

4 .OOOE+OO 

16 

9 . 463E+00 

-4 .OOOE+OO 

17 

0 

1 .437E+01 

18 

-1 .437E+01 

-4 . 545E-07 

19 

-1 .054E+01 

4 .OOOE+OO 

20 

-1 .437E+01 

4 . 545E-07 

21 

-1 .054E+01 

-4 .OOOE+OO 

